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Combining the covariant coalescence model and a blast-wave-like analytical parametrization for 
(anti-)nucleon phase-space freezeout configuration, we explore light (anti-)nucleus production in 
central Au-|-Au collisions at ^/Jnn = 200 GeV. Using the nucleon freezeout configuration (denoted 
by FOl) determined from the measured spectra of protons (p), deutrons (d) and ^He, we find the 
predicted yield of ^He is significantly smaller than the experimental data. We show this disagreement 
can be removed by using a nucleon freezeout configuration (denoted by F02) in which the nucleons 
are assumed to freeze out earlier than those in FOl to effectively consider the effect of large binding 
energy value of “^He. Assuming the binding energy effect also exists for the production of ®Li, ®Li, 

®Li and ®Li due to their similar binding energy values as ^He, we find the yields of these heavier 
(anti-)nuclei can be enhanced by a factor of about one order, implying that although the stable 
(anti-)®Li nucleus is unlikely to be observed, the unstable (anti-)®Li nucleus could be produced in 
observable abundance in Au-|-Au collisions at ^sjviv = 200 GeV where it may be identified through 
the p-'^He (p-'*He) invariant mass spectrum. The future experimental measurement on (anti-)®Li 
would be very useful to understand the production mechanism of heavier antimatter. 

PACS numbers: 25.75.-q, 25.75.Dw 


I. INTRODUCTION 

The quest for antimatter has become one of funda¬ 
mental issues in contemporary physics, astronomy and 
cosmology since the discovery of the positron (the anti¬ 
electron) in cosmic radiation [l| which corresponds to the 
negative energy states of electrons predicted by Dirac 0 . 
Based on very general principles of relativistic quantum 
field theory, it is believed that each particle has its corre¬ 
sponding antiparticle of the same mass (but the opposite 
charge) and any physical system has an antimatter ana¬ 
log with an identical mass. Indeed, following the observa¬ 
tion of antiprotons (p) Q and antineutrons (h) , more 

complex antimatter nuclei such as antideutrons (d) ii, 
antihelium-3 (^He) Q and antitritons (^H) 0 have been 
observed. In terrestrial laboratories, the antihydrogen 
atoms have also been produced H and can even survive 
for a long time in confinement [^. Recently, STAR col¬ 
laboration at RHIC reported the discovery of strange an¬ 
timatter nucleus, the antihypertriton (^H) , and the 

heavier antimatter nucleus antihelium-4 (^He or a) [l^ in 
Au-|-Au collisions. The ALICE collaboration at LHC also 
claimed the observation of "'^He in Pb-|-Pb collisions M- 
A recent review on antimatter production can be found 
in Ref. (T3| . 

The study of antimatter nuclei production in heavy- 
ion collisions is of critical importance for a number of 
fundamental problems in physics, astronomy and cos¬ 
mology. For example, the precision measurement of the 
mass difference between nuclei and anti-nuclei can test 
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the fundamental CPT theorem for systems bound by the 
strong interaction . The measured production rate of 
light anti-nuclei in heavy-ion collisions provides a point 
of reference for possible future observations in cosmic ra¬ 
diation for the motivation of hunting for antimatter and 
dark matter in the Universe [l6l - [T^ . The antimatter nu¬ 
clei production provides the possibility to test the in¬ 
teractions between antimatter and antimatter In 

addition, the production of light anti-nuclei in heavy-ion 
collisions can be used to extract the freezeout informa¬ 
tion of antinucleons in these collisions, which is useful to 
infer the properties of a new state of matter, i.e., quark- 
gluon plasma (QGP) possibly formed in these collisions 
as well as to understand how the QGP expands, cools 
and hadronizes, providing a new window (compared to 
the electromagnetic and hadronic probes) for exploring 
the dynamics of ultrarelativistic heavy-ion collisions. 

The heaviest antimatter nucleus observed so far is ^He, 
and it will remain the heaviest stable antimatter nucleus 
observed for the foreseeable future [l^ barring some dra¬ 
matic discoveries in space detectors due to some special 
production mechanism plj ] or a new breakthrough in ac¬ 
celerator technology. This is because the (anti-)nucleus 
production rate in these heavy-ion collisions is found 
to reduce by a factor of about 10^ for each additional 
(anti-)nucleon added to the (anti-)nucleus according to 
the measured yields of p (p), d (d), ^He (^He) and "^He 
(^He) [l^, and thus the yield of the next heavier sta¬ 
ble antimatter nucleus, antilithium-6 (®Li), is expected 
to be down by a factor of about 10® compared to "^He 
assuming ®Li production rate follows the same exponen¬ 
tial reduction law and is beyond the reach of current ac¬ 
celerator technology. However, a careful observation on 
the experimental yields of p (p), d (d), ®He (®He) and 
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■^He (■^He) in central Au+Au collisions at s/snn = 200 
GeV from STAR [l^ indicates that while the yields of p 
(p), d (d) and ^He (^He) follow an exponential reduction 
rate very well, the yield of "‘He (‘‘He) displays a signif¬ 
icant enhancement (excess) compared to the exponen¬ 
tial reduction rate. The coalescence model calculations 
also significantly underestimate the yield of ^He (^He) 
although they can successfully describe the yields of p 
(p), d (d) and ^He (^He) [l^. It is thus of great interest 
to understand the physics behind this enhancement for 
the yield of ^He ("‘He), which would be critically impor¬ 
tant for the future searching for heavier antimatter nuclei 
such as antilithium-5 (^Li) and ®Li in heavy-ion collisions. 
In Ref. [HI, the possibility of direct production of anti¬ 
matter nuclei out of the highly correlated vacuum has 
been discussed, which provides a potentially more copi¬ 
ous production mechanism for heavier antimatter nuclei 
in heavy-ion collisions. 

In the present work, we propose that the enhancement 
of the ^He (‘‘He) yield could be due to its large binding 
energy which leads to relatively earlier formation for ‘‘He 
(^He) than for d (d) and ^He (^He) in the heavy-ion col¬ 
lisions. Assuming the similar binding energy effects also 
exist for the production of ®Li, ^Li, ®Li and ®Li , we find 
the predicted yields of these heavier (anti-)nuclei can be 
enhanced significantly, implying that although the stable 
(anti-)®Li nucleus is unlikely to be observed, the unsta¬ 
ble (anti-)®Li nucleus could be produced in observable 
abundance in Au-|-Au collisions at RHIC. 


II. THE THEORETICAL MODELS 

Understanding particle production in heavy-ion col¬ 
lision at different energy regions is among the funda¬ 
mental questions in nuclear and particle physics. The¬ 
oretically, the microscopic coalescence model [HI - [H| and 
the macroscopic thermal model (26l - [30| provide two im¬ 
portant approaches to describe the light cluster produc¬ 
tion in heavy-ion collisions. In particular, these two ap¬ 
proaches have been successfully applied recently to de¬ 
scribe the production of light (anti-)nuclei in ultrarel- 
ativistic heavy-ion collisions at RHIC and LHC ener¬ 
gies [^, [29l - [H| . In the present work, the theoretical for¬ 
mulism for the description of (anti-)nuclei production in 
heavy-ion collisions is based on the covariant coalescence 
model [15 together with the (anti-)nucleon phase-space 
freezeout configuration described by a blast-wave-like an¬ 
alytical parametrization [s^ which has been shown to be 
very successful to describe the hadron phase-space freeze¬ 
out configuration in ultrarelativistic heavy-ion collisions. 


A. (Anti-)Nucleon phase-space freezeout 
configuration 

One basic ingredient of the coalescence model is the 
emission source function, i.e., the phase-space freezeout 


configuration, of the constituent particles. In principle, 
the phase-space freezeout configuration can be obtained 
dynamically from transport model simulations for heavy- 
ion collisions (see, e.g.. Refs. [40l - l43| l. In the present 
work, for simplicity, we describe the (anti-)nucleon phase- 
space freezeout configuration using a fireball-like model 
through a blast-wave-like analytical parametrization [s^ ■ 
We assume that particles are emitted from a freezeout 
hypersurface where the particles are in local thermal 
equilibrium described by Lorentz invariant one-particle 
distribution function f{x,p) given by [15 

f{x,p) = gh~^[e^p{{p^u^,-VL)/kT)±l]-^ 

= g(27r)"‘‘[exp(p'"u^/fcr)/^± 1]"‘, (1) 


where the reduced Planck constant ^ ^ is set to be 1, 

g is spin degeneracy factor, H is the chemical potential, 
^ = exp{VL/kT) is the fugacity which is directly related to 
particle number density, is the four-velocity of a fluid 
element in the fireball and T is the corresponding local 
temperature. For the phase-space freezeout configura¬ 
tion, instead of using the four coordinates it, x, y, z), it is 
convenient to use the cylindrical coordinates (r, r, ips, v) 
(see, e.g.. Ref. [15) where r = is the longi¬ 

tudinal proper time, ry = ^lii(|~) is the longitudinal 
space-time rapidity, r is the transverse radius, and (j)s 
is the spatial azimuthal angle. Similarly, the four mo¬ 
mentum {E,px,Py,Pz) is transformed to {mT,PT,4>p^y) 

where pt = \JPx + Py i® transverse momentum, 

y = ^ ln( ) is rapidity, rriT = yjm? -|- p^ is trans¬ 
verse mass, and (j)p is the azimuthal angle in momen¬ 
tum space. The four coordinate and momentum can thus 
be expressed as x^ = (rcoshry,rcos(()s,r sin 0 s,rsinh? 7 ) 
and = {rriT coshy,pT cos (l)p,pT sin (pp,mT sinhy), re¬ 
spectively. 

The one-particle invariant momentum distribution can 
be obtained as 


E 


dp^ 


d^A^ 

prdprdpdf/fp 


= / d^crpp'"/(a:,p). 


( 2 ) 


with dp^ = UprdpTdpd^p. For the particle production 
at midrapidity in heavy-ion collisions that we are con¬ 
sidering in this work, we adopt the longitudinal boost 
invariance assumption [d^ . By setting the longitudinal 
flow velocity vl = zjt, the longitudinal flow rapidity 
gfiow = 5 lii[(l + J^l )/(1 — *^l)] will be identical to the 
space-time rapidity rj = ^\n[{t + z)/{t — z)], and thus the 
four-velocity can be expressed as 


= cosh p{r, (j)s ) (cosh? 7 , tanh p(r, cjig ) cos 4>b , 

tanhp(r, ijis) sin (/)b, sinhp), (3) 


where p is the transverse rapidity of a fluid element in 
the fireball. The above expression can be obtained by a 
longitudinal boost with velocity tanh 77 multiplied by a 
transverse boost with velocity tanh p [39l| . If we fix the 
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freezeout hypersurface by choosing a constant proper 
time T, i.e., is independent of the transverse coordi¬ 
nates, then the covariant normal vector can be expressed 
as 


d^cr^ = (cosh(77), 0,0, — sinh(?7))rrdrd?7d(/)s- (4) 
Therefore, one can obtain 


and (j)h = 4>a 1 and thus the invariant distribution function 
can be expressed as 


f{x,p) = 


9 

(27r)3 


exp ((tot cosh p cosh(?7) — 


Pt sinhpcos((;ip - (j)s))/kT)/^ ± 1 


( 12 ) 


p^Ufj^ = rriT coshp cosh{r] — y) — Pt sinkp cos{(j)p — (j)b), 
p^d^CT^ = ttot cosh(p — p)d77rdrd^s, (5) 

where ipb is azimuthal direction of the transverse flow [3^ . 
Considering that the freezeout can happen in some time 
interval, a Gaussian distribution for the freezeout proper 
time is introduced as follows 


J(t) 



2(At)2 


which satisfies 


( 6 ) 


/ 


J(T)dT = 1, 


/ 


rJ(r)dr = Tq, 


where tq is the mean value of r and Ar is the dispersion 
of the T distribution function. Therefore, the momentum 
distribution can be obtained as 

(fN f 

—-— = / tot cosh(ry-y)/(x,p)J(r) 

PT^PT^y^(pp J 

rdrdprdrd^s. (7) 

Following Ref. [1^ , we parameterize the transverse ra¬ 
pidity of a fluid element in the fireball as 


p = Pol'il + ecos(2((ib)], (8) 


where po is the isotropic part of the transverse flow, e 
is the anisotropic part, ((if, is azimuthal direction of the 
transverse flow which is not identical to spatial azimuthal 
angle and f is the “normalized elliptical radius” 


[rcos((is]^ [rsin(()s]2 

Rl ^ Rl 

with 

R 

ian(t>a = tan ((if,, 

where = i?o(l + S 2 ) is the minor axis of the ellipse, 
Ry — Rq{1 — S 2 ) is the major axis, and S 2 is the geometric 
anisotropy. Therefore, the transverse rapidity can also be 
written as 



(9) 

( 10 ) 


P — POi 


[rcos^s]^ [r sirups 


[Rl 


i ?2 

lly 


-[1-b ecos2((if,]. (11) 


For midrapidity region (y = 0) in central heavy-ion 
collisions that we are considering here, one has S 2 = e = 0 


with p = pot/Rq. One can thus use formula Q and 
(HU to calculate transverse momentum distribution of 
midrapidity particles in central heavy-ion collision. 


B. Covariant coalescence model 

In this work, we calculate light (anti-)nucleus produc¬ 
tion in ultrarelativistic heavy-ion collisions using the co¬ 
variant coalescence formulism [^. In the coalescence 
model, the probability for producing a nucleus is deter¬ 
mined by the overlap of its Wigner phase-space density 
with the nucleon phase-space distribution at freezeout. 
We consider that M nucleons are combined to form one 
nucleus and the total multiplicity of the nucleus can be 
obtained as 

„ M 

= 9c (Y[dNi)p^{xi,...,XM;Pl:---,PM) 

f ( ^ \ 

= 9c f nrfT*J(r*)pfd^(Ti^-^/(xi,pj) j X 

pY{xi,...,xm;Pi,-,Pm), (13) 

where p'^{xi, ...,xm',Pi, ■■■,Pm) is the Wigner density 
function which gives the coalescence probability, pc is 
the coalescence factor [U. By inserting S function to 
conserve momentum, the invariant differential transverse 
momentum distribution of the nucleus becomes 

M 

pY{ xi, ...,xm]Pi, ...,Pm)S^{P - ^Pi). (14) 

The above formula is Lorentz invariant [ 3 ^ and the 
Wigner function is a Lorentz scalar. 

For the Wigner function pY{xi, ...,xm',Pi, ■■■,Pm ), fol¬ 
lowing Ref. [iO, 113, instead of calculating it directly 
using four dimensional coordinators Xi and four dimen¬ 
sional momenta pi of the constituent nucleons, we cal¬ 
culate it in the rest frame of the nucleus. To do so, a 
Lorentz transformation is performed to obtain the space- 
time and energy-momentum coordinates of each nucleon 
in the rest frame of the nucleus. To determine the spa¬ 
tial coordinates of the nucleons at equal time in the rest 
frame of the nucleus, i.e., ri, r 2 ,..., rjv/, the nucleons that 
freeze out earlier are allowed to propagate freely with 
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constant velocity given by the ratio of their momentum 
and energies in the rest frame of the nucleus, until the 
time when the last nucleons in the nucleus freezes out. 
Furthermore, in order to calculate the Wigner function, 
Jacobi coordinate is adopted by a transformation of the 
coordinate as follows l47l | 


/ R \ 


/ ri \ 


qi 


r2 



= Jm 


, (15) 


\ qAf_i / V tm / 


where 


R = 


1 m, 


is the center-of-mass position vec¬ 


tor of the nucleus and ^ — Ti+i) is 

V ^j=i 

the relative coordinate vector. Correspondingly, in the 
momentum space, one has 


f ^ \ 


/ Pi \ 


ki 


P2 



II 


(16) 

V kM-1 / 


V Pm / 



where P is the total momentum of the nucleus and is 
the relative momentum vector. The determinant of the 
Jacobi matrix is \Jm\ = and one then has the 

following identity 

M M-1 

n d^a^.d^p, = d^Rd^P Y[ (17) 

i=l i=l 


Furthermore, we assume the harmonic wave function 
for all the light (anti-)nuclei in the rest frame except the 
(anti-)deutrons for which we use the well-known Hulthen 
wave function (see, e.g.. Refs. [43,|4l|). The Wigner func¬ 
tion of the nucleus can then be obtained as 

pY (xi, ...,xm;Pi, 

= ■ -,(11^1-1, ki,- ■ ■,kM-i) 

M-l 

= 8^”^ exp [ - ^ {qf/af + cr^fc^)], (18) 

i=l 

with af = {miw)~^ where the harmonic oscillator fre¬ 
quency Lu is related to the root-mean-square (rms) radius 
of the nucleus as follows 


4) 


3 1/w 

2=1 


M 

E 

,, 2=1 


rui 



(19) 


Therefore, af can be determined by (rf^) 


of nil = 1^2 = 


2M 


^ IJl \ 
3{M-1) V M/- 


In the case 
= nriM = ’’m, one can obtain 


TABLE 1: Parameters of the blast-wave-like analytical 
parametrization for (anti-)nucleon phase-space freezeout con¬ 
figuration. 



T(MeV) 

Po 

Ro (fm) 

To (fm/c) At (fm/c) 



FOl 

111.6 

0.98 

15.6 

10.55 3.5 

10.45 

7.84 

F02 

111.6 

0.98 

12.3 

8.3 3.5 

21.4 

16.04 


III. RESULT AND DISCUSSION 

A. (Anti-)Nucleon freezeout configuration from 
light (anti-)nuclei production 


We focus on the midrapidity light (anti-)nuclei pro¬ 
duction in central AuJ-Au collisions at ^Jsnn = 200 
GeV in this work. In this case, there are totally six pa¬ 
rameters in the blast-wave-like analytical parametriza¬ 
tion for (anti-)nucleon phase-space freezeout configura¬ 
tion, namely, the kinetic freeze-out temperature T, the 
transverse rapidity po, the longitudinal mean proper time 
To, the time dispersion At, the transverse size at freeze¬ 
out i? 0 : and the fugacity of particle 

For proton phase-space freezeout configuration, we ob¬ 
tain the local temperature T = 111.6 MeV, the transverse 
rapidity po = 0.978, and a constraint on the combination 
of the proton fugacity ^p, tq. At and Rq, by fitting the 
measured spectrum of protons in Au-f Au collisions at 
y/sNN = 200 GeV for 0-5% centrality [i^. To extract 
the values of ^p, tq, At and Rq, we further fit the mea¬ 
sured spectra of deuterons and ^He [i^ simultaneously 
using the results from the coalescence model (see the Sub¬ 
section UTIB] for the details), which leads to Rq = 15.6 
fm. To = 10.55 fm/c. At = 3.5 fm/c and ^p = 10.45. For 
antiprotons, we assume they have the same phase-space 
freezeout configuration as protons except the fugacity is 
reduced to = 7.84 to describe the measured yield ratio 
p/p = 0.75 |48l |. Table U summarizes the parameters of 
the blast-wave-like analytical parametrization for (anti- 
)nucleon phase-space freezeout configuration (denoted as 
FOl). It should be pointed out that we have neglected 
the difference between protons and neutrons (antiprotons 
and antineutrons) for the phase-space freezeout config¬ 
uration due to the small isospin chemical potential at 
freezeout in Au-f Au collisions at ^snn = 200 GeV [1^ . 
Based on the freezeout configuration of (anti-)nucleons, 
one can then predict the production of light (anti-)nuclei 
using the coalescence model. 


B. The production of light (anti-)nuclei 

We use the coalescence model described above to cal¬ 
culate the production of light (anti-)nuclei. In the coales¬ 
cence model, the statistical factor pc is quite important 
and it is given by pc = [l^ with j and N being, 

respectively, the spin and the nucleon number of the nu- 
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TABLE II: Statistical factor Qc, root-mean-square radii 
J'rms [ 13 , HH and binding energy Ei, [s^l of light (anti-)nuclei 



d(d) 

‘He (‘He) 

"‘He (^He) 

®Li (‘Li) 

‘Li (®n) 

5c 

(2xl-t-l) 

-22- 

(2xi + l) 

-23- 

(2x0+1) 

24 

(2x§ + l) 

-23- 

(2X1-I-1) 

2^ 

rrms (fm) 

1.96 

1.76 

1.45 

2.5 

2.5 

Eb (MeV) 

2.224 

7.718 

28.296 

26.330 

31.994 


cleus. The spins of d, ^He, '^He, ®Li and ®Li are 1, 1/2, 
0, 3/2 and 1, respectively. Furthermore, the rms radius 
^rms of the light nucleus is also important since it de¬ 
termines the harmonic oscillator frequency parameter uj 
in the Wigner function of the nucleus. The rrms of d, 
^He, ^He, ®Li and ®Li are taken to be 1.96 frn, 1.76 fm, 
1.45 fm, 2.5 fm and 2.5 fm, respectively [13, [5l|. Here 
the rrnis = 2.5 fm for ®Li is estimated based on the work 
in Ref. [5ll| . For the antinuclei, we assume they have 
the same ground state properties as their corresponding 
nuclei. Table m summarizes the statistical factors, rms 
radii as well as the binding energies [11] of different light 
(anti-)nuclei. It should be mentioned that while d (d), 
^He (^He), ^He (^He) and ®Li (®Li) are stable, ®Li (®Li) is 
unstable against the proton (antiproton) decay with half- 
life of about 370 x 10“^"* s (i.e.. Ill fm/c) [1^ and thus 
it may be identified through the p-^He (p-^He) invariant 
mass spectrum in heavy-ion collisions. 

Figure [T] shows the predicted midrapidity transverse 
momentum distributions of p, d, ^He, '^He, ^Li and ®Li 
together with the experimental data of p from PHENIX 
collaboration [i^ and the data of p, d, ^He and '^He from 
STAR collaboration [ill4i,[ll in central Au-|-Au colli¬ 
sions at = 200 GeV. It is seen that the coalescence 

model predictions with the freezeout configuration FOl 
are in very good agreement with the measured transverse 
momentum distributions of p, d and ^He as expected but 
significantly underestimate the measured yield of ^He by 
a factor of about 6. The similar feature was also observed 
in the calculations in Ref. [^ . 

From Table m one can see that "‘He has a specially 
larger binding energy value compared to d or ^He, and 
thus it is more tightly bound and could be formed in 
relatively earlier stage in heavy-ion collisions compared 
to d or '^He. Physically, the light nuclei can be formed 
in principle in the whole dynamical process of heavy-ion 
collisions, but they are usually destroyed immediately af¬ 
ter their formation due to the violating collisions in the 
high temperature environment. However, for the light 
nuclei with large binding energy values such as "‘He, the 
survival probability in relatively earlier stage in heavy- 
ion collisions is expected to enhance compared to the 
loosely bound d and ^He. In principle, these effects can 
be studied using transport model simulations with dy- 
namic light cluster production in heavy-ion collisions |55l | 
although this is highly nontrivial and beyond the scope 
of this work. In the present work, to effectively mimic 
this binding energy effect, we assume the volume (time) 



FIG. 1: Transverse momentum distributions of light nu¬ 
clei at midrapidity (y=0) in central Au-|-Au collisions at 
-^snn = 200 GeV predicted by coalescence model with FOl 
(solid lines) and F02 (dashed lines). The experiment data of 
protons is taken from the PHENIX measurement whereas 
those of light nuclei are from the STAR measurement [T3, [l^j . 
The data point of protons from STAR measurement has been 
scaled by a factor of 0.6 to correct the weak decay effects [s^. 


of the freezeout hypersurface for nucleons coalesced into 
"‘He is smaller (shorter) than that of d or ^He. For sim¬ 
plicity, we reduce R and t by a factor 1.27 to fit the 
measured yield of "‘He, and this leads to i? = 12.3 fm, 
T = 8.3 fm/c. At = 3.5 fm/c and = 21.4, which is 
denoted as the phase-space freezeout configuration F02 
and is summarized in Table HI Since the binding energy 
values of ^Li and ®Li are also large and comparable with 
that of "‘He as shown in Table El it is thus expected 
that the nucleons coalesced into ®Li and ®Li should have 
similar phase-space freezeout configuration as those coa¬ 
lesced into ■‘He. In Fig. [H we also include the predicted 
transverse momentum distributions of ^He, ®Li and ®Li 
with the freezeout configuration F02. It is seen that us¬ 
ing the freezeout configuration F02 significantly enhance 
the yields of ®Li and ®Li by a factor of about 9 and 16, 
respectively, compared to the case using FOl. 

Figure [2] shows the differential invariant yields 
{(PN/{2'KpTdpTdy)) of (anti-)nuclei evaluated at the 
transverse momentum pt/\B\ — 0.875 GeV/c as a func¬ 
tion of baryon number B. One can see that the coa¬ 
lescence model with FOl reproduces the measured dif¬ 
ferential invariant yields of p (p), d (d) and ’‘He (‘He) 
very well but significantly underestimates the measured 
value of "‘He ("‘He), as already observed in Fig. [TJ The 
dashed lines in Fig. [2] are obtained by fitting the differ¬ 
ential invariant yields of p (p), d (d) and ‘He (‘He) by 
using an exponential function It is seen that the 

differential invariant yields of p (p), d (d) and ‘He (‘He) 
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FIG. 2: The differential invariant yields d?N/{2'KpTdpTdy) of 
(anti-)nucleus at the transverse momentum pr/\B\ = 0.875 
GeV/c as a function of baryon number B in central Au+An 
collisions at ^/snn = 200 GeV. The solid (open) triangles 
represent the coalescence model predictions with FOl (F02). 
The data point of protons from STAR measurement has been 
scaled by a factor of 0.6 to correct the weak decay effects [s^ . 


follow the exponential function very well, depicting the 
same exponential reduction rate of the differential invari¬ 
ant yields with the increased atomic mass number for p 
(p), d (d) and ^He (^He). For ^He ("^He), however, the 
measured differential invariant yields significantly devi¬ 
ate from the exponential function although the coales¬ 
cence model prediction with FOl still follows the same 
exponential reduction rate. For ®Li (®Li) and ®Li (®Li), 
the predicted differential invariant yields with FOl de¬ 
viate the exponential reduction rate by an enhancement 
factor of about 2.25 (3.57) and 1.87 (3.31) whereas those 
with F02 display a much stronger enhancement by a fac¬ 
tor of about 16.6 (26.3) and 23.3 (41.2), respectively, in¬ 
dicating a very strong binding energy effect. 

Table Hill lists the pr-integrated yield in the midrapid¬ 
ity region (—0.5 < y < 0.5), i.e., dN/dy at y = 0, of 
light (anti-)nuclei. One can easily obtain the midrapid¬ 
ity yield ratios d/p (d/p)= 4.65 x 10“^ (3.47 x 10“^) 
and ^He/d (3lfo/d)= 1.99 x lO-^ (1.50 x lO-^). I^par- 
ticular, if we take the midrapidity yield of '^He ("^He) as 
the value predicted by the coalescence model with F02, 
we find the midrapidity yield ratios of ^Li/'*He (^Li/"'^He) 
and ®Li/'*He (®Li/^He) increase, respectively, from 1.32 x 
10-3 (0.99x 10-3) and 1.67x 10-® (0.94x 10-®) with FOl 
to 11.7 X 10-3 (8.74 X 10-3) and 26.4 x 10-® (14.9 x 10-®) 
with F02. These results indicate that the binding energy 
effects can enhance the midrapidity yields of ^Li (^Li) and 
^Li (^Li) by a factor of about 8.77 (8.82) and 15.9 (15.8), 
respectively. 


TABLE III: pT-integrated yield in the mid-rapidity (—0.5 < 
y < 0.5) of light (anti-)nuclei. 



P 

d 

3He 

"He 

®Li 

®Li 

FOl 

16.1 

7.49E-02 

1.49E-04 

1.54E-07 

1.22E-09 

1.53E-12 

F02 

16.1 

- 

- 

9.18E-07 

1.07E-8 

2.43E-11 


P 

d 


"Hi 

®Ll 


FOl 

12.1 

4.21E-02 

6.29E-05 

4.88E-08 

2.88E-10 

2.73E-13 

F02 

12.1 

- 

- 

2.91E-07 

2.54E-09 

4.32E-12 


IV. CONCLUSION 

Based on the covariant coalescence model with a 
blast-wave-like analytical parametrization for the (anti- 
)nucleon phase-space freezeout configuration, we have ex¬ 
tracted (anti-)nucleon freezeout information in central 
Au-I-Au collisions at y/SMN = 200 GeV by htting the 
measured spectra of protons, deuterons and 3He. We 
have found that the covariant coalescence model with 
the obtained (anti-)nucleon phase-space freezeout config¬ 
uration significantly underestimates the measured yield 
of ^He (^He). We have shown the predicted ^He yield 
can be enhanced to the measured value by using a nu¬ 
cleon freezeout configuration in which the nucleons are 
assumed to freeze out earlier than those coalesced into 
deuterons and 3He to effectively consider the large bind¬ 
ing energy value of “^He. The similar conclusion has been 
obtained for ^He. 

Assuming the similar binding energy effect also exists 
for the production of heavier (anti-)3Li and (anti-)3Li due 
to their comparable binding energy values with ^He, we 
have predicted the spectra and yields of (anti-)3Li and 
(anti-)®Li in central Au-|-Au collisions at ^Jsnn = 200 
GeV. Our results indicate that the binding energy ef¬ 
fect can significantly enhance the yields of (anti-)3Li and 
(anti-)®Li. In particular, the midrapidity yield ratios 
^Li/'^He (^Li/'^He) and ^Li/'^He (^Li/'^He) increase, re¬ 
spectively, from 1.32 x 10-3 (0.99 x 10-3) and 1.67 x 10-® 
(0.94x 10-3) without binding energy effects to 11.7x 10-3 
(8.74 X 10-3) and 26.4 x 10-® (14.9 x 10-®) with bind¬ 
ing energy effects. Our results imply that although the 
stable (anti-)®Li nucleus is unlikely to be observed, the 
unstable (anti-)®Li nucleus could be produced in observ¬ 
able abundance in ultrarelativistic heavy-ion collisions at 
RHIG where it may be identified through the p-"^He (p- 
^He) invariant mass spectrum. 

Our present study suggests that the future experimen¬ 
tal measurement on the production of (anti-)®Li in cen¬ 
tral Au-f An collisions at y/SNN = 200 GeV would be 
extremely useful to test the binding energy effect on the 
light (anti-)nuclei production, and thus to understand 
the production mechanism of heavier antimatter nuclei 
in ultrarelativistic heavy-ion collisions, especially the ob¬ 
served enhancement for the yield of ^He (“^He) compared 
to those of p (p), d (d) and 3He (3He). Any deviation of 
the measured (anti-)®Li yield in central Au-f An collisions 
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at ^/snn = 200 GeV from the coalescence model predic¬ 
tion with or without considering the binding energy ef¬ 
fect may indicate the exist of new excitation mechanism, 
e.g., the direct production of nuclei out of the highly cor¬ 
related vacuum. 
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